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Abstract 


An advective mixed-layer ocean model has been 
developed by eliminating the assumption of horizontal 
homogeneity in an already existing mixed-layer model, and 
then superimposing a mean and anomalous wind driven current 
field* This model is based on the principle of conserva- 
tion of heat and mechanical energy and utilizes a box grid 
for the advective part of the calculation. Three phases of 
experimentation have been conducted. In the first phase, the 
model’s ability to account for climatological sea-surface 
temperature (SST) variations in the cooling and heating 
seasons was evaluated. These experiments showed that the 
inclusion of advection results in a large improvement of the 
model’s accuracy and also illustrated the relative importance 
of advection and mixed-layer depth variations. The second 
phase of experimentation consisted of a series of sensitivity 
tests in which the effect of hypothetical anomalous winds was 
evaluated. These tests showed that sustained highly anomal- 
ous winds are capable of generating large-scale SST anomalies 
In the third phase of experimentation, a thirty-day synoptic 
calculation with the model was conducted. For the case 
studied, the accuracy of the predictions was improved by the 
inclusion of advection, although non-advective effects appear 
to have dominated. 
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1. Introduction 


It is v/ell known that the upper layers of the 
oceans undergo substantial temperature variations and that 
these variations are primarily brought about by solar and 
long-wave radiation, heat exchange with the atmosphere and 
heat transfer within the ocean itself. Despite numerous 
studies cf oceanic heat transfer, there are still conflict- 
ing opinions regarding the role of temperature advection 
and the need for its inclusion in numerical models of the 
upper ocean. . 

Hamias (1959) was among the first to suggest that 
anomalous temperature advection in the ocean could be a 
major, cause of large scale sea-surface temperature (SST) 
anomalies. Utilizing the principle that anomalous winds pro 
duce an anomalous drag on the surface water and force an 
anomalous Ekman drift, he computed surface water displace- 
ment vectors by applying Ekman’ s empirical formula 


V = 0 .0127 
W ' / sin 0 


( 1 ) 


to the seasonal anomalous geostrophic wind distribution, where 


Vo is the speed of the surface water current, W is the surface 
wind speed and 0 is the latitude. He then superimposed these 
displacement vectors on the normal sea-surface temperature 
field in order to compute simple advect.ive changes. These 
advectively computed SSI anomalies were then compared with 
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the observed anomalies. Although only slight agreement 
between the two fields could be seen, Namias was able to 
conclude that advection does contribute to the formation 
of SST anomalies but that in this case, other factors must 
have been dominant. He also suggested that anomalous winds 
could produce anomalous upwelling, downwelling and air-sea 
heat exchange. However, no attempt was made to determine 
the relative importance of each. 

Following this early work of Namias, Eber (1961) 
used monthly rather than seasonal wind anomalies and the SST 
fields for the months being studied instead of long-term 
normals. For the same data period as Hamias (1959), he found 
a much closer agreement between the advectively computed SST 
anomalies and the observed anomalies for the winter case, 
while for the- fall case he found only a slight agreement. 

This led him to conclude that advection was the dominant fac- 
tor for the winter studied although other factors must have 
dominated during the fall. He then suggested that seasonal 
variations of the thermal structure in the wind-mixed layer 
might be significant but did not attempt to evaluate its effect. 

Subsequent empirical studies by Hamias (1969, 1970, 
1971, 1972, 1973, 1974) have added further credence to the 
idea that advection within the oceans can be a major cause of 
SST variations, and therefore the generation of large-scale 


SST anomalies. In a study of large-scale variations in sea- 
surface temperatures in the North Pacific, Namias (1970) 
illustrated how water masses could be transported around 
the North Pacific gyre and hovr anomalous winds or currents 
could affect this transport. In 1953 the zonal westerlies 
were stronger than normal and the torque on the Pacific gyre 
was increased, while in 1963 the zonal westerlies and sub- 
tropical easterlies were much weaker than normal and the 
torque on the Pacific gyre was greatly decreased. (This 
torque was actually reversed during the winter of 1963.) 
Namias linked these anomalous winds to the large scale SST 
anomalies of these years and also presented correlations of 
the zonal wind and sea- surface temperatures v/ith season. A 
strong negative correlation was noted for the winter case. 

In his evaluation of the physical causes of the 
winter 1971-1972 SST anomalies, Namias (1972) presented diag 
nostic computations of the anomalous air-sea heat exchange 
which occurred during this period. These values revealed 
that anomalous heat exchange could not have produced the ob- 
served anomalous SST variations and that other factors such 
as advecti on and mixed-layer depth changes must have been 
dominant. However, as noted by Namias (1972), there was in- 
sufficient data to evaluate the effects of changes, in the 
mixed-layer depth* 


As a result of these empirical studies by tfamias 
and those of Bjerknes (1966, 1969), the role of the oceans 
as a lower boundary in atmospheric prediction models has re- 
ceived greater emphasis and it has become increasingly im- 
portant to predict SSI variations accurately. Atmospheric 
response to SSI anomalies is one of the major problems of 
meteorology today and it is generally accepted that accurate 
long-range weather forecasting is not possible unless SST 
variations can also be predicted* 

Initial attempts at extended and long-range numer- 
ical weather prediction utilized models in which the SSI 
field was specified cliraatologically and did not vary with 
time. A series of experiments in which SSI anomalies were 
superimposed on the climatoligical SSI fields, was conducted 
independently by Eowntree (1972), Spar (1973), and Houghton, 
et al« (1974), to determine the influence of large-scale 
SSI anomalies in such atmospheric general circulation models* 
Ihese experiments showed that the model atmosphere is indeed 
sensitive to large, persistent SSI anomalies, and that re- 
mote as well as local effects can be significant. 

later experiments by Spar and Atlas (1975) and Spar 
Atlas and Kuo (1975) , in which synoptic rather than climat- 
ological SSI fields were specified in an atmospheric model, 
indicated that the use of specified synoptic SSI data or the 


daily updating of such fields does not necessarily result 
in an improvement of the atmospheric predictions. This is 
partially due to the loss of predictability in such models 
(Sommerville et al, 1974, Druyan et al, 1975, Seidman 1975) 
and the poor quality of synoptic SST data (Saur 1963, Kao et 
al. 1972, Jastrow and Halem 1973, Smith et al. 1974, Rao 
1974) . . ■ 

Despite these somewhat negative results, the pre- 
diction of SST variations remains as one of the most im- 
portant problems in long-range weather forecasting, and 
such predictions would no doubt benefit the maritime and 
fishing industries as well. Although empirical methods 
have been devised (Naraias 1968, 1969, 1972, 1973, Hammond 
1974) to utilise, SST predictions for long-iange forecasting, 
the major thrust of recent research has been to develop 
predictive ocean models which could he coupled to atmos- 
pheric general circulation models . The predictions for 
both the atmosphere and oceans could then be carried out 
simultaneously, allowing the predicted changes in each med- 
ium, to affect the other. It would thus be possible to study 
the complex interactions which result from anomalies in 
both the atmosphere and ocean. 

The purpose of this particular study is to devel- 
op an advective mixed-layer ocean model hy incorporating the 
effects of horizontal temperature advection into an already 
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existing mixed-layer model, and thpn to study the role of 
temperature advection in the prediction of sea-surface tem- 
perature variations. While this model was designed for the 
purpose of eventually coupling it to an atmospheric model, 

l 

no attempt has been made to do the coupling in this study. 
Instead this study will serve to demonstrate the effect of 
advection by mean ocean currents and indicate how anomalies 
in the atmospheric circulation can result in SST variations. 
This latter point should be extremely important for predic- 
tive runs as well as for simulation studies with coupled 
atmosphere-ocean models. 

Specific objectives of this research are to deter- 
mine (1) if advective effects are significant within the 
mixed-layer ocean model, (2) if the inclusion of advection 
results in an improvement in the accuracy of the SST predic- 
tions, (3) for which regions (and times) advective effects 
are most important, (4*) the relative importance of advection 
and mixed-layer deepening in predicting sea-surface temper- 
ature changes, and (5) the effect of anomalous wind- generated 
ocean currents on the development and maintenance of SST 
anomalies within the model. In order to carry out these 
objectives, the advective scheme must include the advective 
effects of anomalous wind-generated currents as well as mean 
ocean currents. 
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After a brief review of ocean modeling, the de- 
velopment of the advective mixed-layer ocean model and 
results of experimentation with this model will he presen- 
ted. Three phases of experimentation were conducted in 
order to carry out the above objectives. These consisted 
of? (1) climatological simulation, in which the model's 
ability to predict climatological SST variations was eval- 
uated, (2) sensitivity tests, in which the effect of anom- 
alous winds was determined, and (3) a synoptic calculation, 
in which an attempt was made to predict the observed thirty 
day SST change from January 1 to January 31 * 1974 - 


2, Review of Ocean Models 


2.1 Background 

As stated by Miller (1973) , there are two primary 
approaches to the development of ocean models suitable for 
coupling to atmospheric models. She first of these is the 
development of baroclinic models of the entire (three- 
dimensional) ocean which are based on the principles of 
conservation of mass, momentum, energy, heat, and salt. 

In these models, the finite difference approximations to the 
Kavier-Stokes system of equations (Bryan 1975) are simpli- 
fied by assuming that (1) density differences in the ocean 
are negligible in comparison with the mean density 4 , and 
therefore the mean density can be substituted for the actual 
density everywhere but in the buoyancy term (the Boussinesq. 
approximation), (2) the vertical scale in the ocean is much 
smaller than the horizontal scale, and therefore the vertical 
acceleration terms are negligible in comparison with the 
buoyancy term, (3) molecular viscosity is negligible, (4) 
the depth of the ocean is small in comparison with the earth' 
radius, and (5) the angular velocity of the ocean's motion 
is much smaller than the earth's angular velocity. Various 
other simplifying assumptions are included depending on the 
particular model. 


In recent years, a number of such general ocean 
circulation models (Bryan 1969, 1975, Cox 1975, Takano 1975) 
have been developed. While models of this type are desirable 
for the study of oceanic processes and indeed necessary for 
the long-term integrations of climatic studies, they in 
general do not adequately resolve the upper layer of the 
oceans and also have very long response times when coupled 
with atmospheric models. Ihe second approach to ocean model- 
ing is to model the active upper layer of the oceans only. 
Since these models have short response times when coupled 
vdth atmospheric models, they are considered to be preferable 
for short and medium-range predictions. 

Iwo basic approaches to modeling the upper ocean 
exist. Ihe first approach is to account for SSI changes 
wholly in terms of radiation, heat exchange with the atmos- 
phere,' and heat transfer within the oceans. Ihese models 
are either two or three-dimensional in character and do not 
account for temperature changes due to mixed-layer deepening 
or the formation of new mixed layers (i.e. The mixed-layer 
depth does not vary with time). Ihe second approach is to 
develop a model which is capable of predicting mixed-layer 
depth changes. Models of this type are generally one-dimen- 
sional in nature and therefore do not adequately account 
for horizontal heat transfer processes within the oceans. 


2.2 Models With Constant Mixed-layer Depth 


All of the models of this type, to he discussed 
herein, are based primarily on the application of the 
thermodynamic energy equation to the upper layer of the 
ocean and utilize Ekman’s empirical formula (1) to relate 
drift currents in the upper ocean to the low level atmos- 
pheric winds. As mentioned earlier, Namias (1959) was 
among the first to apply Kkman’s formula (1) in an empir- 
ical study of SST anomalies. His model assumed that 
anomalous SSI changes resulted entirely from anomalous 


horizontal advection and therefore, 

3T 


3t — ”"^a 


VSL, , 
m 


( 2 ) 


where Ts is the ambient sea-surface temperature, Tm repre- 
sents the climatological mean sea-surface temperature, V 
is the horizontal gradient operator ( VTm is the horizontal 
mean temperature gradient), and ^a refers t:o the anomalous 
current velocity derived from the use of (1). Since this 
model explicitly neglects all heat sources and sinks, it is 
not suitable for realistic prediction of the SST field. 

Arthur (1966) suggested that an additional advec- 
tion term, representing the adveetion of anomalous isotherms 
by the mean monthly surface currents be included in this 


model. Jacob (1967) incorporated this second advection 
term and an anomalous heating term (derived from empirical 
formulas), into Namias ' model and performed Ms calcula- 
tions over a constant mixed-layer depth of 50 meters 
throughout his grid. Ihese experiments indicated that the 
inclusion of these terms resulted in a significant improve- 
ment in the SSI predictions for a two week to one month 
period. However Jacobs 1 model was not able to account for 
all of the large SST anomalies which were observed. 

Adem (1970) developed a model for the entire 
Northern Hemisphere (with the exception of the Indian Ocean) 
wMch was based on the application of the principle of con- 
servation of mass as well as the thermodynamic energy 
equation. His model included mean and anomalous advection, 
horizontal turbulent diffusion, heat sources and sinks and 
vertical exchange between the upper ocean and the region be- 
low the thermo cline. He assumed that the mixed-layer depth 
was constant at either 50 or 100 meters and performed his 
calculations for both values. His experiments demonstrated 
the usefulness of applying mean and wind drift currents to 
the SST field in order to account for advection as well as 
the relative importance of horizontal temperature advection 
in this model. 

Clark (1972) performed experiments with a model 
which included heat source and sink terms and anomalous 


advection only. He computed SSO? anomalies that would result 
from the inclusion of only advection, evaporative heat ex- 
change, or total surface heat exchange individually, {Dhese 
experiments showed that the adVectively computed SSH? anomal- 
ies are closely related to atmospheric pressure anomalies 
and that anomalous advection can be an important factor in 
bringing about anomalous 3ST variations. However, as indica- 
ted earlier by Eber (1961) and Namias (1972) , and as stated 
by Clark: (1972) the omission of mixed-layer depth changes 
in these calculations is a serious deficiency and therefore 
such changes should be Included in any predictive model or 
specification procedure for SST variations. 

2.3 The Development of Mixed-layer Models 

Although predictive models of the wind-mixed 
layer are relatively recent in origin, attempts to model 
and/or describe mixed-layer depth changes have been going 
on since Ekman’s (1905) study of the influence of the earth's 
rotation on ocean currents, Ekman assumed that the water 
was homogeneous and applied a constant value of the eddy 
viscosity coefficient to obtain his famous "spiral'’ calcu- 
lation. This "spiral" concept therefore did not contain a 
well mixed upper layer bounded by a distinct gradient below. 


However, Ekman realized that his theory was an oversimpli- 
fication and stated that the coefficient of eddy viscosity 
cannot be constant when the density of the water is not 
uniform (as is the case in the real ocean). 

Munk and Anderson (1948) developed a density 
stratified model of the upper ocean utilizing eddy coeffic- 
ients of viscosity and conductivity which varied with depth. 
They obtained a steady state (non-time dependent) solution 
for the distribution of temperature and velocity with depth 
which accounted for the upper mixed layer and the thermo- 
cline below. However, their computations yielded thermo- 
cline depths which were in general less than half of the 
observed depths and they were only able to apply their theory 
where evaporation was small, the heat flux was positive, and 
advection was negligible. 

Other steady state models of the thermocline have 
been developed by Kitaigorodski in I960 and Kraus and Eooth 
in 1961 (see Kraus and Turner, 1967) . Kitaigorodski computed 
the depth of the mixed layer from the balance between the’ 
mean work of the wind stress and the work required to mix 
incoming heat downward, while Kraus and Rooth showed how the 
absorption of heat through a finite depth in conjunction with 
a net heat loss at the surface (through evaporation, conduc- 
tion, and long-wave radiation) resulted in convection within 
the mixed layer. In recent years several time-dependent, 



one-dimensional mixed-layer ocean models have been developed 
for the purpose of short and medium-range SST prediction. 

These models -which are primarily based on the conservation 
of heat and mechanical energy, as well as suitable boundary 
conditions, provide reasonable predictions for the changes 
in mixed-layer depth and sea-surface temperature, hut be- 
cause of the assumption of a horizontally homogeneous ocean 
in these models and their one-dimensional nature, they neglect 
or do not adequately taKe account of important dynamical 
effects such : as advection of heat hy horizontal ocean cur- 
rents and upwelling and downwelling. 

She first realistic time-dependent model of the 
seasonal thermocline was developed by Kraus and (Turner (1967) . 
They assumed that heat fluxes into the mixed layer and mass 
entrained at the bottom of the layer are instantaneously 
mixed uniformly throughout the layer, and thereby considered 
the mixed layer to be vertically homogeneous. As a result 
of this assumption, it was not necessary for them to utilize 
eddy coefficients in their formulation. Kor the case of con- 
stant mixing by light winds (approximately equal to 5 meters 
per second) , they were able to obtain quantitative as well 
as qualitative predictions for variations in sea-surface 
temperature and mixed-layer depth. It was found that their 
model’s predictions agree reasonably well with oceanic ob- 
servations and with earlier laboratory experiments (Turner 
and Kraus 1967) • : 
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In a subsequent laboratory experiment, Kato and 
Phillips (1969} studied the development of a mixed layer in 
an annular tank, They applied a constant stress to the 
surface of an initially quiescent fluid with a uniform den- 
sity gradient, and observed the growth of a turbulent layer 
by entrainment of the underlying fluid. It was found that 
the rate of increase of the potential energy of the stratified 
fluid is proportional to the rate of dissipation of kinetic 
energy in the turbulent layer, and that the thickness of 
this turbulent layer increased as the time to the one-third 
power, Kato and Phillips reasoned that these results should 
be directly applicable to the real ocean. 

Penman (1973) developed a time-dependent, one- 
dimensional model of the upper ocean, which is essentially a 
generalization of the Kraus and Turner model and is consis- 
tent with the results of Kato and Phillips. This model 
assumes that the ocean is a stably stratified, incompressible 
fluid obeying the Boussinesq approximation, and it ignores 
wave-induced dynamical effects. The ocean is considered to 
be horizontally homogeneous except that indirect effects can 
result from the non-zero curl of the wind stress. 

figure 1. illustrates the fundamental parameters 
and boundary conditions of the Penman model. Prom this figure, 
it can be seen that the vertically homogeneous mixed layer is 
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bounded on the bottom by a temperature d is continuity t below 
which a climatological temperature gradient is specified. 

£he -vertical temperature profile may be specified there- 
fore in terms of three variables: !P S . (the temperature of 

the mixed layer) , T-h (the temperature immediately below 
the mixed layer), and h (the mixed-layer depth). 

Prom the conservation equations for thermal and 
mechanical energy, Denman derived the following set of three 
first order ordinary differential equations for the predic- 
tion of these variables: 

9 

ci!E 

\ I - (G-D)+h'(B+H +H )+R{h-Y~ 1 +y- 1 e~ yh ] , (3) 

h 

,, 2 [G-D+Ry" 1 (l-e~ Th ) ] -h [B+H +H C +R (l+e“ Yh ) ] 

H [ wtffi] = — — ' e -' > - - , (4) 

h(E s -T_ h ) 



YR e -1fh - (w+ah/at) 


-h 


(5) 


In the above equations, G represents the production of tur- 
bulent energy from the wind stress at the sea surface, D is 
the total energy dissipation within the mixed layer, y is 
an average extinction coefficient for solar radiation, w is 
an imposed vertical velocity at the bottom of the layer, and 
H is a Heaviside step function which is equal to sero if 
w + dh/dt * 0 or one if w + dh/dt >0* 2Ke terms 3, He, Hs 


< ■ - i J .... .1 1 1 . 
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and R represent the surface fluxes of the net long-wave 
radiation# latent heat, sensible heat, and short-wave solar 
radiation, all divided hy pc where P and c are the den- 

P P 

sity and specific heat of sea water respectively, i.e. 

* 

(H, B, H e> H s ) = r|- <R*, B*, He*. H s *) 

P 

where the * denotes that these fluxes are at the sea surface. 

Che set of coupled equations ( 3-5) is similar to 
the predictive equations of Kraus and Turner’s (1967) model. 
However, Kraus and Turner solved their model analytically 
on a time scale of months whereas this model is solved numer- 
ically on a time scale of days. Because of this, the boundary 
conditions for the Denman model have been formulated somewhat 
more precisely. Specifically, the Denman model includes a 
specified temperature gradient below the temperature discon- 
tinuity and allows the penetration of solar radiation into 
this region. 

Obviously the predictions of the above system of 
equations depend very strongly on the value of the Heaviside 
step function H. If H = 1, then the model is in a wind dom- 
inated regime and the entrainment mixing term (w + dh/dt) is 
included in ( 4 ) . Under these conditions, the mixed layer is 
deepening and cold water from below the temperature discon- 
tinuity will be entrained into the mixed layer, due to the 


t 
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work done "by turbulence against the buoyancy. The tempera- 
ture of the mixed layer would tend to decrease due to this 
effect. 

When H = 0, the model is in a heat-dominated re- 
gime and there is no entrainment mixing term in ( 4 ) . Since 
the absorption of solar radiation is greatest near the sur- 
face of the ocean and all of the available turbulent energy 
from the wind stress must he used to redistribute this heat 
uniformly throughout the mixed layer, none of the energy is 
available to deepen the mixed layer further. Therefore, in 
this case, the mixed-layer depth will remain unchanged or 
under conditions of rapid heating and low winds, a new shall- 
ower mixed layer of warmer water will form. This new mixed 
layer would be superimposed on the old temperature profile 
and there would thus be more than one temperature discontin- 
uity in the profile. 

Henman performed a number of experiments with this 
model for both the wind-dominated and heat-dominated regimes. 
For the wind-dominated case, he found that (1) doubling the 
turbulent energy available for mixing over a two day period, 
produced a significant increase (almost 30^) in the mixed- 
layer depth, (2) decreasing the stratification in the model 
to one half of its original value resulted in a smaller hut 
still important effect on mixed-layer depth, (5) typical 


summer heating rates (solar radiation of 400 cal. cm 4 ^ day~^ 
coupled with back radiation of -80 cal. cm, day ) did not 
significantly affect the mixed-layer depth, and (4) large 
evaporative heat losses (typical of winter conditions) 
coupled with high winds cause strong convective mixing which 
does have an important effect on the mixed layer, For the 
heat dominated case, he found that the predictions of the 
model are sensitive to the value of the extinction coeffic- 

f 

ient and that decreasing this extinction coefficient to one 
half of its original value increases the mixed-layer depth 
signif icantly (70$ in the case studied) . He also found that 
if the extinction term was not retained in the region below 
the mixed layer, that the mixed-layer depth was overestimated 
by 15$. 

Denman and Miyake (1973) utilised the Denman model 
to predict SST and mixed-layer depth changes at Ocean Station 
Papa during a twelve day period from 13 to 24 June, 1970. 
Oceanographic observations, described in their paper, indi- 
cate that during the summer, sensible and latent heat fluxes, 
horizontal temperature advection, and vertical advection from 
belov/ the mixed layer may all be neglected at this location 
(50 K, 145 V/) in the Pacific Ocean. Under these conditions, 
the formation of a shallow mixed layer of warm water is caused 
primarily by intense solar radiation and low wind speeds. 
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Deepening of the layer is induced primarily by the increased 
winds accompanying atmospheric storms. For this period the 
Denman model's predictions of SST agreed well with observa- 
tions and the model's profiles did simulate the time-depend- 
ent behavior of the mixed-layer. 

Two other one-dimensional models have been developed 
to predict SST variations and to explain the time-dependent 
behavior of the mixed-layer, Pollard, Rhines and Thompson 
(1973) developed a model, which like Denman' s, describes the 
response of the upper ocean to an imposed wind stress and 
heat flux. In the Denman model, part of the available tur- 
bulent energy from the wind stress induces a slow erosion of 
the stably stratified, quiescent layer below the mixed layer. 
However, in the Pollard, Rhines and Thompson model, available 
turbulent energy is used to drive inertial motions within the 
mixed layer. These motions allow for a more rapid erosion of 
the stably stratified layer and .therefore a more rapid mixed- 
layer deepening. Hiiler (1973) has developed a model which 
includes both of these erosion mechanisms. He states that 
both of these processes are important at different times. 
However the necessity for including inertial motions for SST 
prediction is not well established. 

As previously mentioned, there has been a major 
effort in recent years to develop coupled atmosphere-ocean 
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models, and for short and extended-range prediction, it is 
desirable to utilize a model of the upper ocean for this 
purpose. At the Goddard Institute for Space Studies (G-ISS), 
the one dimensional mixed-layer ocean model developed by 
Denman has been adapted for global use in conjunction with 
the GISS atmospheric model by Miller (1974). This GISS ocean 
model has a horizontal grid spacing of 4° in latitude and 5° 
in longitude. Despite the fact that this is a global grid, 
the one-dimensional nature of the Denman model has not 
changed and the predictions at each point are carried out 
independently. Some modifications to the original Denman 
model have, been made by Miller in order to improve the SSI 
and mixed-layer depth predictions, These will be discussed 
briefly in the next sections. 

As a result of the initial experimentation with 
this GISS ocean model and for the purposes listed in the in- 
troduction, an advective mixed-layer ocean model has been 
developed by incorporating advective effects into this model. 
In order to account.for the horizontal advection of heat at 
each gridpoint, the assumption of a horizontally homogen- 
eous temperature field for the entire ocean has been elim- 
inated (thus allowing for horizontal temperature gradients) 
and a mean and anomalous wind driven current field has been 
superimposed on the grid. By talcing this approach we still 


allow the modified Denman model to account for the effects 
of solar and hack radiation, sensible and latent heat ex- 
change with the atmosphere, and mixed-layer depth changes 
at each point, while the superimposed current field serves 
to couple the gridpoints through the advection of heat and 


mass . 


3.» Description of the Advective Mixed-Layer Model 
3.1 Basic theory 

With, the exception of the advection scheme, the 
model to "be described herein is almost identical to the 
Denman model and therefore our approach will he very sim- 
ilar, As in the Denman model, we assume that the ocean 
is a stably stratified, incompressible, fluid obeying the 
Boussinesq approximation and ignore wave induced dynamical 
effects, molecular heat fluxes and the viscous generation 
of heat. We do not however, assume that the entire ocean 
is horizontally homogeneous in all properties at this time, 
{Thus our derivation will include horizontal advection terms 
for heat and mass. Horizontal turbulent diffusion will be 
assumed to be negligible when compared to the other heat ex- 
change processes (Clark 1972) and will not be included. It 
can be argued (Bathen 1971) that lateral diffusion is im- 
portant in certain regions. Although this is probably true, 
these effects will not be considered here. 

Figure 1 depicts the main parameters and processes 
of the Denman model. The GISS ocean model profile is only 
slightly different from this and is depicted in four forms 
in Figure 2. 


z-joom 


Schematic illustration of the vertical 
temperature profile of the G-ISS ocean 
model at four different stages of develop 
ment . (a)' Initial "profile. ( b ) Typical 
profile at end of cooling season. ( c) 
Typical profile at beginning of heating 
regime, (d) Typical profile during the 
heating season. ^ sl >^ S 2 and -s4>~s3 >T s2 
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When the model is initialized, the initial temperature pro- 
file at a gridpoint will he similar to that shown .in Figure 
2a, This structure consists of a vertically homogeneous 
mixed layer hounded hy a temperature discontinuity, below 
which two temperature gradients are specified. The upper 
temperature gradient 
thermocline while the lower one max corresponds to 

the temperature gradient below 100 meters. This structure 
is somewhat more realistic than Denman's for modeling an 
entire cooling season. 

During the cooling season the mixed layer deepens, 
and by the end of this season the mixed-layer depth is gen- 
erally between 100 and 200 meters in mid latitudes (Miller 
1974) . The solid line in Figure 2b represents a typical pro- 
file at the end of the cooling season (Miller 1974) . The 
deepening of the mixed layer is arbitrarily cut off at a 
maximum depth of 200 meters. If this point is reached, then 
the mixed layer continues to cool until the temperature jump 
at the bottom goes to aero or a heating regime begins. When 
the heating regime begins, a new (shallower) mixed layer forms 
at the surface and this is superimposed on the preceding tem- 
perature profile. In this case, there is a new isothermal 
layer and temperature discontinuity (as shown by the dashed 
line in Figure 2b) . The temperature gradient immediately be- 
low this is initially set equal to the remaining portion of 
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corresponds to the seasonal 
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the old isothermal layer (i*e. — - 0, initially)* The 

initial heating regime profile that would result from this 
is represented by the solid line in Figure 2c. As the heat- 
ing regime continues a shallower mixed layer forms and this 
new layer is superimposed on the preceding profile. This 
is illustrated by the dashed line in Figure 2e. At this 
stage, the profile is geometrically adjusted to assure that 
there are equal areas of warming and cooling on each side 
of the upper temperature discontinuity. This eliminates 

this temperature discontinuity and yields a new temperature 
3T 

gradient -g— k , as illustrated by the dotted line in Figure 
2c. The complete profile which would result from this pro- 
cess is depicted in Figure 2d, 

The adjective model’s profile is exactly the same 
as that described above except that horizontal ocean cur- 
rents are included in the mixed layer. These currents are 
averaged within the mixed layer, and for simplicity it Is 
assumed that the current speed goes to zero at the bottom 
of; the mixed layer. 
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3.2 Model Equations 


Poliowing Denman (1973) the thermodynamic energy 
equation for the ocean may he written as 


dT/dt = Q T /pc p 


(6) 


where IP is the temperature and is the heat source, due 
to absorption of solar radiation. Denman (1973) makes the 

Y2 

assumption that Qrj, = YR* e j where z is the negative depth 
in the ocean. If we substitute this expression into the 
time averaged, Eulerian, turbulent form of ( 6 ) we obtain 


3T 

at 


+ V * 7T + 

' 3z nr ' 


yz 


Pc, 


( 7 ) 


(1) (2) (3) 

where 'J is the horizontal current velocity and the primes 
in this and the following equations denote perturbation 
quantities, I.e. deviations from local, climatological mean 
values. lerm 1 is the local change of temperature, term 2 
represents the horizontal temperature advection by ocean 
currents, and terra 3 is the local divergence of turbulent 
heat flux. In this equation, the vertical temperature ad- 

flirt 

vection term has been eliminated because of the absence 

of vertical temperature gradients within the mixed layer, 
and horizontal eddy diffusion has been neglected. 
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Denman (1973) employs Phillips (1966) turbulent 
kinetic energy equation. 


JL(°i = _siwT\8£_X r w < (El + si)] -W'-P ’g _ 

2 } 3Z 9Z lw l p Q 2 ; J p Q 


( 1 ) 


( 2 ) 


(3) 


(4) 


where c represents the eddy kinetic energy, S 
the mean horizontal current in the mixed layer, 
pressure, 5 is the gravitational acceleration, 


e r (8) 

represents 
P denotes 
and e re- 


presents the rate of dissipation of turbulent energy. This 
equation is used without modification in the advective model 
on the assumption that the horizontal advection of kinetic 


energy is negligible. 

In equation 8, term 1 represents the time rate of 
change of kinetic energy of the turbulent motion, which 


according to Denman (1973), is negligibly small. Term 2 
represents the rate of production of turbulent energy by the 
turbulent Reynolds stresses, which act on the mean current 


shear. Term 3 is the local divergence of the vertical trans- 
port of turbulent mechanical energy. (It is this term which 
serves to instantly mix the turbulent energy uniformly through- 
out the layer . ) Term 4 is the rate of turbulent energy loss, 
due to work done against the lower density gradient. The 
elimination of term 1 results in a steady state relationship, 


expressing the balance between energy source and sink terms. 
Thus the steady state mechanical energy equation for the 
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turbulent motion is 

2 

-s T w r ||-4-[w'(|-+V)J = “gw'®' + s (9) 

It is now necessary to consider the application 

of the principles of conservation of heat and mechanical 

energy (equations 7 and 9) to the mixed layer , the interface 

below the mixed layer, and the region helow this interface. 

Although Figures 1 and 2 depict an artificial temperature 

discontinuity at Z - -h it is more reasonable to consider 

the interface between the mixed layer and lov/er layer as 

having some thickness Ah , This is illustrated in Figure 

3, which is taken directly from Denman’s (1973) paper. 

Denman (1973) evaluates the mixing entrainment at 

the bottom of the mixed layer ( w'T ' j _^) by integrating the 

thermodynamic energy equation across the interface 

-h+Ah -h+ h z 

f [f£+V ■ VT + w |£+ JL (^V)] - / - dZ (10) 

~h 8t 32 82 dz -h pG p 

(Note that in (id) the vertical advection term has been in- 
cluded since this region is below the mixed layer and verti- 
cal temperature gradients do exist.) The integration above 
results in 


<f£+v VT) Ah + (T s -T_ h ) (w+^) t w'T’ 


3h, 


-h+Ah 


_ YR*P 


-yh 


P Q cp 
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If we neglect horizontal variations in h and allow Ah-K> 
as depicted in Figures 2 and 3 , then following Denman (1973), 

“ -H(W+||} (T s -T_ h ) (12) 


where H is the Heaviside step function (defined previously)* 
At the top of the mixed layer, the turbulent heat 
flux is eq.ua! to the net heat transfer across the ocean sur- 
face. ffltius, 


wT — -F 
o 


at the upper "boundary, where 


F = H s + H e ..+ B. 


(13) 


(14) 


Within the mixed layer itself, following Denman 
(1973), integration of the thermodynamic energy equation (.7) 
with the aid of 12 and 13 results in 


.ST 


h(- 


+ V • +H(w+fe {T -T,) = P + R(l-e" ,n ) , (15) 


-Yh, 
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where V is the vertically .integrated current in the mixed 
layer, i.e , 


V = t / vd'Z 
h -h 


(16) 


Integration of the steady state mechanical energy equation 
( 9 ) within the mixed layer results in 


•/°S T w T |fdz - W (J“r + % ? -> 
-h 8z po 2 


z=Q 


ag w>r 32 + Edz 
-h 


(17) 


where it has been assumed that 


w'(£l+%-) = o . 

p 0 Z=-li 


Prom (17) it follows that 
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If we let 
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and G- = ~ 
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then 


0 


-/ w'T'dZ = G- D 
-h 


Poliowing Denman (1973), if we now integrate 
equation 7 twice with respect to Z and combine it with 
equation 20 we arrive at 

grjT 

+ v • VT = (G-D) + Fh+ R(h -y" 1 ) 4- Y~ 1 Re “ Th 

h 


V/e can eliminate the terms on the left hand side of (21 ) 
hy combining equations 15 and 21. 3}his results in 

H (w + « 2 [G-D-FRy " 1 (l~e“ yh ) ] - h [F+R(l+e" Yh) J. - 


h (T s - T h ) 


( 18 ) 


(19) 


( 20 ) 


( 21 ) 


( 22 ) 
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for tile region below s » -h, while at Z = -h 


~(T~h) = yRe~ Th - 


/ , dh . 8T 
dt az 


-h 


(24) 


Here it has been assumed that no turbulent energy penetrates 
below z=-h and as mentioned earlier no currents exist at 
Z *s -h or below. 

We can now restate the complete set of equations 
for the adveotive mixed-layer model as 


3T 


VT = ~[-(G-D)+h(B+He+Hs')+R(h-Y 1 +y“ 1 e“ ,Yh ) ] 
it 


flrw+ & = 2 [G-D+Ry" 1 (l-e" Th ) ] -h (B+He*Hs+R (1+e 
L dt J h(T s -T h ) 


-yh 


)] 


d m = v R e" Yh - ( w + ~) — 

T -h ^ Ke dt j az 


dt 


(25) 


(26) 


(27) 


~h 


As can be seen, this system of equations is almost 
identical to Denman’s (equations 3-5 in this paper) and dif- 
fers primarily in the inclusion of the horizontal adveotive 
term v • VT in equation 25 . All' of the parameters for the 
above model are readily available from routine meteorological 
and oceanographic data with the exception of G- and D. Denman 
assumes that we may approximate the term G - D as follows 

G - D = -255- . (28) 

P 0 ag • 

In the above equation, m is a constant fraction of the rate 
of transfer of turbulent energy downward from the wind field 
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at 10 meters, and E a the rate of work done "by the wind stress 

"3 

at 10 meters is equal to P a c ^0 U iQ ' where P a is the den- 
sity of the air, C^q is the drag coefficient at 10 meters 
and the mean wind at 10 meters. Miller (1975) has 

found that this formulation results in too much deepening of 
the mixed layer during the cooling season and has empirically 
modified (23) so that for the cooling season only 


G-D = 


iriEa 

p 0 ag 


.4 h~ 
dt 


(29) 


where Q 

A = / T d Z . 

~h 

She addition of this empirical factor results i‘h a more real- 
istic prediction for the change in potential energy of the 
mixed layer and hence the deepening of this layer during the 
cooling season. 

We may solve the system of equations (25 - 27) hy 
first separating equation 25 into advectfve and non-advective 
parts. If we assume that there is horizontal homogeneity in 
the immediate vicinity of each gridpoint hut that the entire 
ocean is not horizontally homogeneous, then as previously 
noted, we may utilize the modified Denman model (without the 
advective terms) to account for the SST changes due to solar 
and hack radiation, air-sea heat exchange, and mixed-layer 
depth changes at each gridpoint. The advective temperature 
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change may then he computed independently at all gridpoints. 
The combination of these two calculations results in the 
final SST prediction for each gridpoint. 

This process may he stated more precisely hy the 
following steps: 

1. ) Solve equations 25-27 at each gridpoint wif! the term 

V * VT omitted. These one-dimensional calculations re- 
sult in independent predictions for T h, and at 
each gridpoint. 

2. ) Solve the advection scheme (to be discussed in the next 

subsection) for all gridpoints. This results in an ad- 
ditional heat flux term for the mixed layer, 

5.) Add the predicted adveotive temperature change to the 

mixed layer calculation at each gridpoint and adjust the 
final prediction of T s , h, and to account for the 

addition or removal of heat from the mixed layer due to 
horizontal advection. 

3-3 The Advection Scheme 

The advection scheme to be described herein con- 
sists of (1) specifying an ocean current field for the entire 
horizontal grid, (2) solving the advective parts of the 
equations for conservation of heat and mass for the entire 
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horizontal grid, and (3) determining the change in tempera- 
ture at each gridpoint due to the net heat and mass trans- 
fers at these points. 

In Section 1, it was stated that the advection 
scheme must include the effects of both mean ocean currents 
and anomalous wind- generated currents, in order to carry out 
the objectives of this study. As a generalization of the 
approach used by Adem (1970) we assume that, 


v V M v AD V MD ; V G 


(30) 


In the above equation the total ocean current V has been 
divided into its mean and anomalous parts, where V M is the 
mean monthly or seasonal ocean current (specified from cli- 
matology, V AD is- the drift current computed from the actual 
(observed or predicted wind), v is the drift current com- 
puted from the mean monthly or seasonal winds, and is an 
anomalous geostrophic current, computed from anomalous changes 
in the sea surface slope. Prom the definition of the ver- 
tically averaged current in the mixed layer (equation 16 ) 

we may define 

- 1 0 ' 

¥ = r / V dZ 

M E M 


V MD h V MD d 2 


? AD = K £ d Z 


(31) 


*= ^rr 
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ana v = ? M + (V AD - v MQ ) + v G (32) 

Thus the total current through the mixed layer consists of 
"both mean and anomalous components* The fields of V M which 
were used in this study were prepared 1 "by the author and are 
described in the next section. As mentioned, this field is 
specified and does not change with time. represents 

the pure wind driven part of V M and also does not change 
with time. V on the other hand, represents the pure wind 
driven current due to varying atmospheric winds and there- 
fore at any instant of time ( V^ D - V MD ) represents the 
anomalous wind driven ocean current. 

Both of these drift currents are computed from 
Ekman theory. According to this theory (Adem 1970 or Neumann 
1968), whenever a wind blows over the ocean surface, pure 
wind drift currents are generated by the transfer of momen- 
tum from the wind stress to the water and by the friction 
between moving layers of water. Surface drift currents in 
the Northern Hemisphere are directed 45 ° to the right of the 
wind direction while each succeedingly lower layer is di- 
rectly at an increasing greater angle and slower speed than 
the layer immediately above it. The level at which the cur- 
rent speed has decreased to e _7r of its surface value and its 
direction is exactly opposite to the surface current direction 
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is referred to as the depth of friction influence. Within 
the layer of frictional influence, above the depth, the 
total current transport vector is directed 90° to the right 
of the wind and its magnitude is equal to 0.225 times the 
magnitude of the surface current. These principles may be 
expressed by the following formulas: 


V DX “ V 0 e“ (ir/D,Z =os< | 


-(ir/D)S . . 7T _TT , 

V DY + V 0 e Sin( 4 D Z) 


(33);. 


where V DX is the component of the pure drift current at a 
depth Z, which is perpendicular to the surface wind direction, 
is the component of the pure drift current at a depth z, 
which is parallel to the surface wind direction, V Q is the 
speed of the surface current, and D is the depth of friction- 
al influence. At the surface Z = 0 and 


TT 


V DX = V 0 COS 4 


V DY ~ V 0 Sin 4 


(34) 


while at the depth of frictional influence Z = D and 


V DX = COS ( l 


V DY - Vq©" 11 sin ( —nr) 


( 35 ) 
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Empirical values for the surface current speed V Q and the 
depth of frictional influence D have "been determined hy 
Ihorade (see Adem 1970). For surface wind speeds greater 
than 6 meters per second, 

v . 0^0127 „ 

/sin <j) 

, C36) 

(which is the same as equation 1) 

and D = 

/sin <j> 

For winds speeds less than or equal to 6 meters per second 

„ 0^ 02 59 ^ 

0 /sin J 

(37) 

and D _ 3.67 SW 

/sin <j> 

From the above formulas we may now state our equa- 
tions for evaluating the pure mean or actual wind drift 
velocity components through a depth h (equal to the mixed- 
layer depth) as 

u_ = C 7 - (in- cos 0 + v TT sin 0) for W > 6 m sec 1 

/sm <j> (38 

v~ = C ( v cos 0 - u TT sin 0) for W > 6 m sec 

D Ssin~J W W 

or • 

Uj- “ C • (u-- cos 9 + vl, sin 0) for W - 6 m sec 1 

/sin "<jT 


v = g (v TT cos 0 ~ in. sin 0) for W - 6 m sec ^ 

D /iliTT W W 
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where and are the west-east and south-north components 
of the mean or actual drift current through the mixed layer, 
u and v v are the mean or actual wind components in the west- 
east and south-north directions, 0 is the angle between the 
wind, and the total wind driven current through the mixed- 
layer, and C is a current magnitude factor. In order to 
determine the parameters 0 and 0, it is necessary to com- 
pute the depth of frictional influence D from equations 36 
or 37 and compare this value to the mixed-layer depth h. 

If h is greater than or equal to h then we assume following 
Ekman that the drift current speeds below the depth D are 
negligibly small and we set C =: 0.225 and 6 - 90°, If h 
is less than D then 0 will lie between 1 and 0.2.25 and 0 
will he between 45° and 90°. The exact value for C and 0 
depends on the ratio h/D. 

Namias, (1959) Eber t (l961) and Jacob (1967) used 
values for C and 0 of 1 and 45° respectively while Clark 
(1972) and Adem (1970) tested their model with values of 
0.225 and 90° as well. They found that their best results 
were obtained with values of C = 1 and 0 - 45° • However 
the author feels that it is not reasonable to consider the 
surface current (as indicated by these values for C and 0 } 
to extend throughout the entire mixed-layer depth. The 
approach taken here has been to allow convergence or diver- 
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gence of anomalous drift currents to result in changes in 
the height of the sea surface at each primary gridpoint. 
Prom the sea surface slope so created, we may compute the 


anomalous geostrophic current V . The components of this 

y 

geostrophic current u and v are given by the following 

y y 

formulas: 


u = — g,— M 
g 2Q sin(j> AY 


v = § M 

g 2£2 smtj) AX 


(40) 


(ALL 

where AZ is the change in the height of the sea surface, 

AZ 

and represent the anomalous sea surface slope, and ^ 
is the earth’s angular velocity. The combination of the 
geostrophic and wind driven components results in a total 
current through the mixed layer, which has a larger magni- 
tude and smaller angle of deviation from the wind than the 
original drift current. This approach is consistent with 
Ehman's "elementary current system" (see Neumann 1968} and 
provides a somewhat sounder theoretical basis for the advec- 
tion scheme. 

The grid system used in this model is a simplified 
version of the "box method" (Kurihara 1967) and has the pro- 
perty of conserving both heat and mass for. the advective 
calculation. Specifically the grid consists of a series of 
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grid boxes, at the center of which are the primary grid-points 
of the GISS model (described earlier) . A section of this grid 
is depicted in Figure 4. Bach individual grid box is horizon- 
tally homogeneous in all properties* However, the entire 
ocean is not assumed to be homogeneous. As previously men- 
tioned, this approach enables us to apply the modified Denman 
model at each of the primary gridpoints depicted by dots in 
Figure 4 to predict T g , h, and T_ h non-advectively. The ocean 
current values are computed at each of the secondary grid- 
points depicted by circles in Figure 4 * Averaging of these 
current values at the four corners of each box results in the 
net zonal and meridional currents for that box* This in turn 
is used to account for the net heat and mass transfer through 
each box. In this manner, heat and mass are exchanged between 
neighboring grid boxes and the temperature advection calcula- 
tion for the entire grid is accomplished. 

If we consider each primary gridpoint to be speci- 
fied by the index i,j (where the i value represents the longi- 
tude, and the j value represents the latitude) , then we may 
specify the four secondary gridpoints surrounding the primary 
gridpoint f s (i+j, j+j) »■ {i+^v j-|) , Ci“|, j-|) and (i-|, jt|), 

and the four sides of the gridbox as (i+|, j), (i, j-|), (i-|>j) 
and (i, - This specification Is depicted in Figure 5* 

The mass transport across each grid box side can be computed 
from the following quantities : 



Fig. 4. A section of the grid for the advective 
cahculation. Dots represent primary 
. gri'dpoints . Circles represent secondary 
gridpoints. _ - 



1 


1 


46 






+ u. ,i . i ) AY . h mean 


2 {n i-bi+i + '"l-i.-H 1 iY 3 h raean 


2 j+£ + ' v i+| ,.j+§^ AX i h mean 


#{V. x -_1 + V. 1 ._i) AX. h mean 

^ ^ 2 r J 2 fj J 3- 


(41) 


where AY^ is the distance between secondary j gridpoints 
(equal to 4° of latitude) , AX^ is the distance between secon- 
dary i gridpoints (equal to 5° of longitude) , and h mean is 
the mean depth of the ocean current across that side of the 
grid box. Summing up these four transports results in the 
net mass transport into the box. If there is net mass con- 
vergence into the box, then AZ (the change of the sea surface 
height at the primary gridpoint) will be positive and the sea 
surface will rise. If there is a net divergence for the box, 
then AZ will be negative and the sea surface will fall. From 
the new values for Z at each primary gridpoint , the anomalous 

geostrophic current y can be computed from equations 40. 

■ g 

These equations can be restated in finite difference form as 
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The thermodynamic energy equation for the advective 
part of the calculation is 


3T ST 3T 

s s s 

_= 


(43) 


where u and v are the west-east and south-north components of 
the total current in the mixed layer. Multiplication of this 
equation by the density p and combination with the continuity 
equation results in the flux form of the advective equation. 
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In order to compute the heat fluxes into the box, 
we must consider the difference in temperature between the 
water inside the box and the water surrounding the box. The 
heat fluxes across each side of the box may be computed from 


(pU Vi + A,j = P U i+i,j ( ‘-‘-2 


T + T. 

S * • S* , | • 
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The local temperature tendency within the "box is then cal- 
culated from the divergence of the heat flux, and the 
temperature at each gridpoint is predicted by the Euler 
method. Stated simply, the predicted temperature at the 
primary gridpoint is equal to the original heat content 
of the box plus the heat flux into the box, divided by 
the new volume of the box. 
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4. Climatological Simulation with the Advective Mixed-layer 
Model 

The initial phase of experimentation with the ad- 
vective mixed-layer model was designed to (1) test the 
model’s ability to predict climatological SST variations, 

(2) determine if the inclusion of advection results in an 
improvement in the predictive skill of the mixed-layer model, 
and (3) determine the relative importance of advection and 
mixed-layer depth changes. Two 3-month prediction experi- 
ments, one in the cooling season and one in the heating sea- 
son, were conducted. We shall refer to these as experiments 
1 and 2 respectively. 

lata for these two experiments consisted of (1) 
initial conditions of the monthly mean SST and mixed-layer 
profiles for each of the primary gridpoints, (2) climatol- 
ogical values for the solar and back radiation, wind speed, 
and sensible and latent heat fluxes, at each primary grid- 
point, (3) climatological ocean current values in component 
form at each of the secondary gridpoints, and (4) monthly 
mean SST fields for verification. All of the SST data for 
these two experiments were derived from the monthly mean 
ocean temperature tabulation developed at the National Center 
for Atmospheric Research (Washington and Thiel 1970), The 




climatological mixed-layer profiles were derived from the 
Mechanical Bathythermograph (MB I) Data Pile of the Nation- 
al Oceanographic Data Center. This data file consists of 
vertical profiles of temperature, with observations at 
five 'meter intervals. The mixed-layer depth was taken to 
be the first point on the profile where the decrease of 
temperature with depth is . 2°C/5 meters or more. 

The mean fluxes were taken directly from a tabu- 
lation of such data by Schutz and Gates (1971). The data 
used by the author to prepare the mean ocean current field 
consisted of (1) climatological observations of surface 
ocean currents, obtained from the Climatological and Ocean- . 
ographic Atlas for Mariners, Volumes 1 and 2 (1959, 1966). 
(2) a limited number of observations of the vertical varia- 
tion of horizontal ocean currents, obtained from the Woods 
Hole Oceanographic Institution and from several standard 
oceanographic texts (Sverdrup et. al, 1942, Defant 1961, 
Stommel 1965, Neumann and Pierson 1966, Stommel and Yoshida 
1972) and (3) a limited number of observations of the mean 
density fields for the Atlantic and Pacific oceans, also 
obtained from the Woods Hole Oceanographic Institution and 
the same standard oceanographic texts mentioned above. These 
sources of data enabled the author to make a subjective ap- 
proximation to ‘the vertical profiles of the horizontal ocean 


current at each gri&point and hence the average horizontal 
current in the mixed layer. Specifically the method followed 
was to (1) interpolate the surface ocean current values to 
the primary gridpoints of the advective model, (2) estimate 
the vertical variation of these currents with depth from the 
Ekman spiral theory, as described in the previous section, 

(3) modify the profiles so obtained using the limited number 
of vertical profile observations and estimates of the geo- 
strophic current shear obtained from the density fields, 

(4) average the vertical profiles to obtain the mean current 
in the mixed layer, (5) divide these currents into their 
west-east and south-north components, and (6) interpolate 
these values to the secondary gridpoints. 

In experiment 1, the prediction was made from Jan- 
uary mean initial conditions and mean fluxes were updated 
daily to drive the model. The mean ocean current values 
for the three month period were specified and remained con- 
stant. Obviously the anomalous current terms (v^ - ) 

and y were equal to zero since the actual and mean winds 
6 ' 

were the same in these experiments. SST predictions from 
three, different versions of the mixed-layer model were eval- 
uated at the end of thirty, sixty, and ninety days. Version 
(A) refers to the original version of the GISS ocean model 
in which no advection is included. Version (B) is the same 
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as (A) except that in this case advection by mean ocean 
currents is included. Version (C) is the same as (B) ex- 
cept that in this case the mixed-layer depth is held con- 
stant with time and therefore no deepening is allowed to 
occur. This last case was included to illustrate the im- 
portance of mixed-layer deepening in the advective model, 
since earlier advective models of the upper ocean (Jacob 
1967, Adem 1970, and Clark 1972) did not account for such 
mixed-layer depth changes. A fourth version, which we shall 
refer to as Version (P), represents the persistence forecast 
and is included for comparison, to determine if any of the 
three forecasts possess skill over persistence. 

SST predictions were made for the entire Northern 
Hemisphere (between 2°N and 70°N) with the exception of the 
Indian Ocean. She results of this experiment and the one to 
follow were analysed in terms of the absolute difference be- 
tween the predicted and observed sea-surface temperature at 
each of the primary gridpoints, for the four different ver- 
sions. These absolute errors were then added for the entire 
North Atlantic and North Pacific (between 2°N and 70°N) and 
then divided by the total number of gridpoints for each of 
these regions to yield the average absolute error for the 
North Atlantic and North Pacific. The average absolute tem- 
perature error for the four versions of the model is presen- 


ted in Figures 6 and 7 for the Forth Atlantic and Forth Pac- 
ific respectively. From these figures, it can he seen that 
the advective effect, as represented by the difference be- 
tween the A and B curves, is relatively large for both oceans 
and that the inclusion of advection has resulted in approxi- 
mately a 50 % reduction in the average absolute error of the 
SSI predictions. The effect of mixed-layer depth changes, 
as represented by the difference between 0 and B is also 
large. In the Forth Atlantic, neglect of mixed layer vari- 
ations (curve C) produces a forecast which is no better than 
that resulting from neglect of advection (curve A), For the 
Forth Pacific, inclusion of advection and neglect of mixed- 
layer depth variations (C) produces a forecast only slightly 
better than that of the non-advective GISS model (A). 

Clearly, both advection and mixed-layer depth variations 
must be included in the ocean mixed-layer model. In order 
to determine the skill of the A, B and C predictions over 
persistence, the average absolute error for these three 
model versions was compared with average absolute error of 
the P forecasts. It is clear from these figures, that al- 
though all three versions possess skill over persistence at 
thirty days, only the version B, which includes both advec- 
tion and mixed-layer depth changes, consistently maintains 
predictive skill over persistence out to ninety days. 


Average Absolute Error (°C). 
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.Big. 6. Average absolute errors in predicted sea- 
surface temperatures (°C) over the North 
■ t - Atlantic, generated by versions (A), (B) , 

(C) and (P) for the cooling season (January 
initial conditions). 
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The significance of the advective effect is also 
illustrated by Figures 8 and 9. Figure 8 presents the error 
field for version A at the end of ninety days while Figure 9 
presents the error field for version B at the end of ninety 
days; A comparison of the figures indicates that virtually 
all ’of the large errors have "been eliminated as a result of 
the advective effect and only a few smaller errors remain. 

Some of these remaining errors are due to errors in the mean' 
ocean current field while others may he due to the simpli- . 
fying assumptions of the mixed layer model. A comparison 
of the absolute errors at each gridpoint for model versions 
A and B shows that for the North Atlantic, the predictions 
were improved at 71$ of the gridpoints, remained the same 
at 10$ and were worsened at 19$ of the gridpoints by the 
inclusion of advection. For the North Pacific, the inclu- 
sion of advection resulted in an improvement of the predic- 
tions at 65$ of the gridpoints, while 22$ were worsened, and 
13$ remained unchanged. Clearly, errors in the mean current 
field are responsible for many of these advective errors. 
However, improvement of the field will probably not occur 
until sufficient observations of the vertical profiles of 
horizontal ocean currents are obtained. From Figures 8 and 
9 it can be seen that the largest advective effects are as- 
sociated with the major currents particularly the G-ulf Stream, 



F£g. -8. Srror field for the 90-day’ sea-surface temperature 
* predictions from January mean initial conditions, 
using model version A. Solid lines are error 
isotherms drawn for an interval of 1 C. 
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Kuroshio, labrador, Canary and California currents. Not 
enough advection of heat has occurred off the west coast 
of Africa or the southern tip of Japan, Y/hile too much ad- 
vection has occurred along the east coast of the United 
States. 

In order to gain further insight into the region- 
al importance of the advective effect, the average absolute 
errors have been computed for the regions of the North At- 
lantic and North Pacific between 34°N and 70°N and between 
2°N and 30°N. She graphs for the regions are presented in 
Figures 10-13. A comparison of Figures 10 and 11 reveals 
that the variability of the climatological SST field, as re- 
presented by the P curve, is considerably larger for the 
northern region of the North Atlantic than for the southern 
region, while the accuracy of the version P predictions are 
similar for both regions. This means that version P possess- 
es somewhat more skill in the northern region than in the 
southern region and that the model is capable of predicting 
both large and small temperature changes. The advective 
effect and the effect due to mixed-layer depth changes are 
about equally large for both regions. However this latter 
effect is actually slightly larger and more important than 
the advective effect at sixty and ninety days for the northern 
region, A comparison of Figures 12 and 13 reveals similar 
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results. However, the effect due to mixed-layer depth 
changes is considerably larger than the advective effect 
for the northern regions of the North Pacific at sixty and 
ninety days. 

In experiment 2, the ninety- day predictions were 
made from June initial conditions and the mean fluxes were 
again updated daily to drive the model. Data for this ex- 
periment was prepared from the same sources and in the same 
manner as in experiment 1. Ihe summer mean current field 
differs from the winter mean current field primarily in the 
fact that the currents are considerably more meridional in 
direction in the North Pacific while being only slightly 
weaker in strength. Also the mixed-layer depths are much 
shallower in the summer experiment and therefore the upper 
level currents, when averaged over the entire mixed layer, 
are stronger, 

• Ihe average absolute errors for versions A, B, C 
and P are presented in figures 14 and 15 for the North 
Atlantic and North Pacific respectively* A comparison of 
these figures with figures 6 and 7 shows that in experiment 2 
the average absolute errors for version P is more than twice 
as large as in experiment 1 at thirty days and more, than three 
times as large at sixty and ninety days. The average absol- 
ute error for versions A, B and C is only slightly larger 
in experiment 2. As a result of this, all three of these 





versions maintain predictive skill over persistence through- 
out the ninety period, for "both oceans. The advective effect 
and the effect due to mixed-layer depth changes are consi- 
derably larger in this experiment. This is due to the shall- 
ower mixed-layer depths and stronger mean currents, as men- 
tioned earlier. 

Figures 16 and 17 illustrate the ninety day error 
fields for versions A and B respectively. From these figures 
it can be seen that the inclusion of advection has decreased 
the predictive errors, although they have not been completely 
eliminated. In the North Atlantic, the predictions at 64 $ 
of the gridpoints were improved by advection, while 14 $ re- 
mained the same, and 22$ were worsened. In the North Pacific, 
the inclusion of advection resulted in an improvement at 63$ 
of the gridpoints, while 21$ were worsened, and 18$ remained 
unchanged. Once again, the largest advective effects appear 
to be associated with the major current systems. Too much 
advection of heat has occurred off the east coast of the 
United States and in the northern North Atlantic while too 
little advection has occurred in the northwestern and north- 
central North Pacific. 

Figures 18-21 present the average absolute errors 
for the four regions, previously defined. These figures show 
that the average absolute error of the P forecast is consi- 


67 





Average Absolute Error (°C) 




Average Absolute Error (°C) 



te E 




derably larger for the northern regions and that this effect 
is more pronounced in the Pacific Ocean. The advective 
effect is also considerably larger in the northern regions 
and the deepening effect is slightly larger than the advec- 
tive effect in all cases except the thirty and sixty-day 
predictions for the northern North Atlantic. 
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5 . Anomalous Wind Sensitivity Tests 

After the completion of the climatological sim- 
ulation experiments, several sensitivity tests were con- 
ducted to evaluate the role of ocean currents generated by 
anomalous winds. Specifically, these tests were designed 
to determine if anomalous advection by such currents could 
be a major factor in the generation of large scale SST anom- 
alies, as suggested by Namias (1972), and to evaluate the 
effects of specific anomalous wind patterns. The climato- 
logical predictions of the preceding experiments were set 
up as control runs upon which various hypothetical anomal- 
ous winds were superimposed. The sensitivity of the model 
to these winds was then measured in terms of the difference 
between the SST predictions of the control run and each test. 

Data for this phase of experimentation consisted 
of exactly the same data fields as described in the preced- 
ing section plus the climatological mean wind fields for 
the cooling season and heating season experiments. The mean 
zonal and meridional wind components were obtained from the 
TJ.S. Navy wind atlas (Chief of Naval Operations 1966). These 
climatological wind components were specified at each of the 
secondary gridpoints and remained constant throughout the 
ninety day predictions. Hypothetical wind components were 


74 


also specified at each of the secondary gridpoints and held 
constant. Y/ind fields used for these tests were obtained 
by (A) reversing the climatological wind directions, (F) 

i 

doubling and (C) halving the wind speeds, (D) doubling and 
(E) halving the zonal wind components, and (?) doubling and 
(G-) halving the meridional wind components. In additional 
tests, the major pressure and wind patterns were shifted 
(H) 12° north, or (I) 12° south, or ( J ) 10° east or (K) 

10° west. Finally (1) the torque on the major anticyclonie 
gyres of the North pacific and North Atlantic was reversed. 

In the latter test, the speed of the westerlies in the 35°N~ 
55°N band was decreased to 0.1 times its mean value, while 
the easterlies in the 20°N-35°N band were decreased to 0.8 
times their mean speed and reversed in direction. This wind 
field is similar to the extreme 1963 winter case in the North 
Pacific as described by Namias (1970). However, since the 
climatological SSI field was used for this sensitivity test, 
no attempt was made to simulate the 1963 winter SST anomalies. 

The results of these sensitivity tests for the cool- 
ing season experiment are presented in Tables 1 and 2. Table 
1 presents the average absolute SST differences (°G) computed 
for the North Atlantic and North Pacific while Table 2 pre- 
sents the corresponding maximum absolute differences. The 
tables show that the SST variations due to anomalous wind 



Table 1. Average absolute differences in computed sea-surface 
temperatures (degrees C) from wind sensitivity tests 
for the cooling season (January initial conditions). 


A.) Reversed wind directions 

30 

60 

90 days 

North Atlantic 

.19 

.37 

.58 

North Pacific 

.28 

.61 

.89 

B.) Doubled wind speeds 




North Atlantic 

7X3 

— -22 

.36 

North Pacific 

.17 

.36 

.59 

C.) Halved wind speeds 




North Atlantic 

.05 

.11 

”713 

North Pacific 

,08 

.18 

.30 

D.) Doubled zonal components 




North Atlantic 

7T3 



.TZ 

North Pacific 

.19 

.37 

.57 

E.) Halved zonal components 




North Atlantic 

7 m — 

.13 

”773 

North Pacific 

.09 

.19 

.29 


i 
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Table 1. continued 


F.) Doubled meridional components 

30 

60 

90 days 

North Atlantic 

.08 

.17 

.28 

North Pacific 

.09 

.19 

.32 

G. ) Halved meridional components 




North Atlantic 

.04 

.09 

.15 

North Pacific 

.04 

.10 

.17 

H.) 12° northward shift 




North Atlantic 

.11 

.22 

”735 

North Pacific 

.18 

.35' 

.54 

X.) 12° southward shift 




North Atlantic 

.10 

.19 

.33 

North Pacific 

.12 

.30 

.51 

J , ) 10° eastward shift 




North Atlantic 


TT5 

“721 

North Pacific 

.09 

.20 

.31 

K.) 10° westward shift 




North Atlantic 

.06 

.12 

.22 

North Pacific 

.09 

.19 

.31 


] 
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Table l. continued 


1.) Retarded torque 

30 

60 

90 days^ 

North Atlantic 

.10 

e 

CM 

* 

.33 

North Pacific 

.16 

.30 

.45 


Table 2. Maximum absolute differences in computed sea-surface 
temperatures (degrees C) from wind sensitivity tests 
for the cooling season (January initial conditions)* 


A.) Reversed wind directions 


North Atlantic 

1.3 

2.6 

4.2 

North Pacific 

2.2 

4.3 

6.9 


P.) Doubled wind speeds 



North Atlantic 

“7S 

— ITT 

371 " 



North Pacific 

1.2 

2.6 . 

4.2 


c.) 

Halved wind speeds 






North Atlantic 

T4 - 

78 — 

1.4 



North Pacific 

.6 

1.3 

2.0 


33.) 

Doubled zonal components 


■ 




North Atlantic 

“78 — 

' 176 — 

3.3 



North Pacific 

1.4 

2.8 

4.0 



■f 

i 
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Table 2. continued 

£.) Halved zonal components 30 60 90 flays 

North Atlantic .4 .8 1.7 

North Pacific .7 1.5 2.3 


P.) Doubled meridional components 


North Atlantic 

.9 

. 

° 

3- j. 

North Pacific 

.9 

2.1 

3. .3 

G.J Halved meridional components 




North Atlantic 

.4 

ITT 

i'.'a' " 

North Pacific 

.4 

1.2 

1.9 

H.) 12° northward shift 




North Atlantic 

“79 

1.8 

573 

North Pacific 

1.2 

2.2 

3.4 

I.) 12° southward shift 




North Atlantic 

75 

1.2 

2-z 

North Pacific 

.9 

2.0 

3.1 

J . ) 10° eastward shift 




North Atlantic 

7 T“ 

^ 

p-g 

North Pacific 

.8 

2.1 

5.2 


Table 2. continued 


K.) 10° westward shift 50 60 90 days 


North Atlantic 

.5 

1.2 

1.8 

North Pacific 

.6 

2.1 

3.2 


L.) Retarded torque 

North Atlantic .6 1.2 T75 

1.1 2.3 3.5 


North Pacific 


generated currents can be quite large. In some regions the 
effect of sustained highly anomalous winds is at least as 
large as the advective effect due to mean ocean currents, 

The largest effect obtained was due to the rever- 
sal of the wind circulation. Although such a reversal of 
winds for the entire North Pacific or North Atlantic is not 
likely to occur^ persistant reversals of the wind over small- 
er regions have been observed (Namias 1970), Since this 
apparently represents the maximum effect of anomalous winds, 
it is useful to look at the regional generation of SST an- 
omalies in the case of reversed winds. The SST difference 
field after thirty days for this test is presented in Figure 
22. Reversal of the winds produces anomalous southerly 
ocean currents oyer the northern half of the North Pacific 
and over the northern and eastern North Atlantic, while an- 
omalous northerly currents are found over the remainder of 
these oceans. From Figure 22, it can be seen that anomalous 
warming and cooling associated with these anomalous currents, 
resulted in several large SST anomalies. Anomalous warming 
is found along the west coast of California where southeast 
winds have replaced northwesterly winds, setting up an anom- 
alous southerly current. The largest SST anomaly however is 
generated in the northwestern North Pacific. A comparison of 
these differences with those computed from climatological 



Fig. 22. sea-surface temperature anomalies (°C) generated after 30- 
days by reversing the wind directions. 






current advection (not reproduced) reveals that this maxi- 
mum SST anomaly occurs in the region where the strongest 
westerly current component has been reversed. Further com- 
parison of the difference field with the climatological 
advection field shows that every large SST anomaly coincides 
with a correspondingly large wind anomaly. Similar results 
were obtained for all of the sensitivity tests although the 
magnitudes of the differences were smaller. 

The same set of sensitivity tests was also carried 
out for the heating season. The summer mean wind field was 
considerably weaker than the winter field. However, the 
wind generated currents, when averaged in the mixed layer, 
were actually slightly stronger in the heating season due to 
the shallower mixed-layer depths. Tables 3 and 4 present the 
average and maximum absolute SST differences for this exper- 
iment. A comparison of these tables with Tables 1 and 2 
reveals that the SST changes due to anomalous wind-generated 
currents is slightly larger in the heating season than in the 
cooling season. However this difference is smaller than the 
difference between the advective effect of mean currents in 


the two seasons 
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Table 3. Average absolute differences in computed sea-surface 
temperatures (degrees C) from wind sensitivity tests 
for the heating season (June initial conditions). 


A.) Reversed wind directions 30 60 90 days 


North Atlantic 

.24 

.47 

.66 

North Pacific 

.31 

.64 

= 95 

B.) Doubled wind speeds 




North Atlantic 

7X2 — 


TT2 

North Pacific 

.20 

.39 

.61 

C./ Halved wind speeds 




North Atlantic 

7ofT~ 

.13 

.21 

North Pacific 

.09 

.19 

.31 

D. ) Doubled zonal components 



: 

North Atlantic 

7T? 

.29 

.43 

North Pacific 

.19 

.41 

.63 

E. ) Halved zonal components 




North Atlantic 

.07 



North Pacific 

.10 

.21 

.32 


Table 3. continued 


F. ) Doubled meridional components 

30 

60 

90 days 

North Atlantic 

.10 

.21 

.33 

North Pacific 

.13 

.26 

.41 

G. ) Halved meridional components 




North Atlantic 

705 

no 

TT6 

North Pacific 

.06 

.13 

.21 

H.) 12° northward shift 




North Atlantic 

715““ 

.30 

.43 

North Pacific 

.19 

.37 

. 56 

I. ) 12° southward shift 




North Atlantic 

.14 

777 

~T40 

North Pacific 

.17 

.35 

.53 


J.) 10° eastward shift 


North Atlantic 

.08 

7T6~ 


North Pacific 

.09 

.21 

.33 


K. ) 10° westward shift 

North Atlantic 709 7TS T2B - 

.10 .22 .35 


North Pacific 


I 


l 


1 


85 


Table 3. continued 


Xi. ) Retarded torque 

50 

60 

90 days 

North Atlantic 

.15 ' 

.30 

.43 

North Pacific 

.20 

.38 

.57 


Table 4. Maximum absolute differences in computed sea-surface 
temperatures (degrees C) from wind sensitivity tests 
for the heating season (June initial conditions). 


A.) Reversed wind directions 


North Atlantic 

1.8 

" 3T5 

““371“ 

North Pacific 

2.4 

4.9 

7.8 

B. ) Doubled wind speeds 




North Atlantic 

1.0 

2.1 

3.3 

North Pacific 

1.3 

3.0 

4.7 

C.) Halved wind speeds 




North Atlantic 

".5 

1.0 

1.6 

North Pacific 

.7 

1.5 

2.4 

D.) Doubled zonal components 




North Atlantic 

ITT — 

— 271 — 

— 572 

North Pacific 

1.4 

2.7 

4.1 


I 


Table 4. continued 


E.) 

Halved zonal components 

•50 

60 

90 days 


North Atlantic 

.6 

1.1 

1.6 


North Pacific 

.7 

1.3 

2.0 

F.) 

Doubled meridional components 




North Atlantic 

1.2 

2.3 

3.5 


North Pacific 

1.2 

2.4 

3.7 

&.) 

Halved meridional components 





North Atlantic 

.6 

1.1 

1.8 


North Pacific 

.6 

1.2 

1.9 

H.) 

12° northward shift 





North Atlantic 

1 . 2 

2.4 

3.3 


North Pacific 

1.5 

2.7 

4.0 

I.) 

12° southward shift 





North A clantic 

.8 

1.7 

2.7 


North Pacific 

.9 

2.1 

3.2 

J.) 

10° eastward shift 





North Atlantic 

.8 

' * ITS — 

“2T3 


North Pacific 

.8 

1.7 

2.8 


Table 4. continued 


K. ) 10° westward shift 

North Atlantic 
North Pacific 

L. ) Retarded torque 

North Atlantic 


30 

60 

90 day : 

.6 

1.2 

2.0 

.7 

1.3 

2.2 


r;o' "" 

271 — 

570 

1.2 

2.5 

3.8 


North Pacific 
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6. Synoptic Calculation 

In the previous experiments it was demonstrated 
that the advective mixed-layer model is capable of account- 
ing for the climatological SST changes in the cooling and 
heating seasons. The sensitivity of the SST predictions 
to anomalous winds was also shown. In the third phase of 
experimentation with the model, a thirty day synoptic cal- 
culation during the cooling season was performed. This 
experiment was conducted in order to determine (1) if the 
inclusion of advection by mean and anomalous currents re- 
sults in an improvement in the accuracy of the model’s 
synoptic predictions, (2) the magnitude of the advective 
effect, and (3) the relative importance of advection by 
mean end anomalous currents. 

Data for this experiment were similar to those 
used in the climatological cooling season experiment, and 
differed primarily in the use of synoptic SST fields and the 
inclusion of anomalous winds. Daily observed SST fields 
were obtained from the U.S. Navy Fleet Numerical V/eather 
Central {FNWC). These are derived from surface ship and 
buoy observations, supplemented by satellite data, and cov- 
ered the North Atlantic and North Pacific from 18°N-66°N. 
Climatological mean winds were obtained from the U.S, Navy 
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wind atlas (Chief of Naval Operations 1966) while synoptic 
winds were obtained froru the Marine Section of the National 
Weather Service (NWS) in New York City* These synoptic winds 
were averaged over the thirty-day period to yield a mean 
wind field for the month. The climatological values for 
solar and back radiation, sensible and latent heat fluxes 
and mixed layer profiles (as described in section 4) we re 
specified at each primary gridpoint, since there were no 
synoptic observations available for these quantities. 

In this experiment, SST predictions were made from 
January 1, 1974 initial conditions and the mean fluxes were 
updated daily to drive the model. The mean ocean current 
values, climatological mean January winds, and mean January 
1974 winds were specified at each of the secondary grid- 
points and remained constant for the thirty-day period. The 
anomalous current terms were then computed at each of the 
secondary gridpoints and added to the mean current value. 

SST prediction^ were made with several different 
versions of the mixed-layer model in order to evaluate the 
relative importance of these current terms. Version (A) 
once again refers to the original version of the GISS ocean 
model, with no advection included, while Version (B) in- 
cludes advection by mean currents only. (There was no Version 
C in this experiment.) Version (D) refers to the case where 
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advection by both mean and anomalous currents is included. 

In Version (E) only advection by anomalous currents is con- 
sidered, and Version (P) is the persistence forecast. In 
addition, three predictions were computed in which only 
the advection, and no other model physics was included. 
Versions (F), (G), and (H) refer to the cases of mean advec- 
tion, mean and anomalous advection, and anomalous advection 
respectively. In these latter computations, the thirty-day 
SST change is due only to the advective effect. The results 
of all of these predictions were analysed in terms of the 
average absolute errors for the North Atlantic and North 
Pacific and are presented in Table 5. 

From this table, it can be seen that the most ac- 
curate prediction was produced by model Version D in which 
the mixed-layer model physics and all current terms have 
been included, A comparison of Versions A, B, D, and E 
reveals that (1) the improvement of accuracy due to the in- 
clusion of mean currents is small, (?) a much larger improve- 
ment results when all current terms are included, and (3) 
the improvement of the predictions due to advection by 
anomalous currents is more important than the improvement 
due to mean currents in this experiment. A comparison of 
Versions F, G, H and P reveals that the anomalous current 
alone yield the most accurate predictions. All model 



Table 5. Average absolute error, °C, in the thirty-day 
SST predictions from January 1, 1974 initial 
conditions. 



Version 

North Atlantic 

North pacific 

Mixed-Layer Physics plus 




No Advection 

A 

*75 

.71 

Mean Advection 

B 

.70 

.68 

Complete Model 

B 

.57 

.59 

Anomalous Advection 

S 

.66 

.64 

No Mixed layer Physics 

- 



Mean Advection 

F 

00 

00 

• 

.95 

Total Advection 

a 

*<v 

OO 

• 

.87 

Anomalous Advection 

H 

.82 

.79 


Persistance 


P 
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versions maintain skill over persistence with the exception 

* 

of Version F for the North Pacific. 

In order to gain further insight into the importance 
of advecticn by anomalous currents, the maps of mean sea 
level pressure for January 1974, predicted anomalous SST 
change due to anomalous advection (Version H), and observed 
anomalous SST changes are presented in Figures 23-25. Figure 
23 shows the isobars of average sea level pressure for January 
1974. In the North Atlantic, the mean circulation was dom- 
inated by an intense Icelandic low. This low was centered 
close to its normal position but had a central pressure 22 
millibars (mb.) deeper than normal. The subtropical high 
pressure belt in the North Atlantic was close to its normal 
position and had a central pressure only 1 mb. higher than 
normal. In the North Pacific, the western lobe of the' 

Aleution low was dominant in January 1974 and was 4 mb. deep- 
er than normal. Meanwhile the eastern North Pacific high 
was centered slightly further east than, normal and had a 
central pressure 4 mb. lower than normal. 

The predicted and observed anomalies of thirty-day 
SST change are presented in Figures 24 and 25. A comparison 
of these fields reveals that there is a remarkable oualita- 
tive agreement between the two fields. Almost all of the 
major areas of anomalous SST change have been predicted. 
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Fig. 23. Mean sea level pressure pattern (mb.) for January 1974 
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Pig. 23. Mean sea level pressure pattern (mb.) for January 1974. 
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Pig. 24. predicted changes of sea-surface temperatures (°C) 

from January 1 to January 31, 1974 using model version H. 
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However , the magnitude of the advective 
anomalies Is in most cases less titan half of the ob- 
served anomalous change. This indicates that although 
anomalous advection was important it v/as not dominant 
in January 1974. 


7. Summary and Conclusion s 

By eliminating the assumption of horizontal 
homogeneity and superimposing a mean and anomalous wind 
driven current field, the effects of horizontal advection 
of heat have been incorporated into an existing, two- 
dimensional, time-dependent ocean mixed-layer model. 
Mixed-layer depths and temperatures in the Worth Atlantic 
and Worth Pacific Oceans have teen computed with the model 
ty means of a numerical scheme which separates the calcu- 
lation into advective and nonadvective parts. This approach 
allows the original mixed-layer model to account for the 
effects of solar and tack radiation, sensible and latent 
heat exchange with the atmosphere, and mixed-layer depth 
changes at each point, while the superimposed current field 
serves to couple the gridpoints through the advection of 
heat and mass. 

r 

Three phases of experimentation have teen conduc- 
ted with this model. In Section 4, the model’s ability to 
account for climatological SST variations was evaluated 
through two experiments - one in the cooling season and 
one in the heating season. These experiments showed that 
(1) the inclusion of advection results in a marked improve- 
ment in the accuracy of the ocean model climatology, (2) 
the largest advective effects are associated with the major 
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current systems, (3) the advective effect is larger in the 
heating season than in the cooling season, and (4) the ad- 
vective effect and the effect due to mixed-layer depth 
variations are of similar magnitude. 

In Section 5, the sensitivity of the model SST 
predictions to various hypothetical anomalous wind distri- 
butions was evaluated. The results of this phase of exper- 
imentation indicate that the ocean model (and presumably the 
real ocean as well) responds sensitively to anomalous winds, 
and that sustained highly anomalous winds are capable of 
generating large-scale SST anomalies. The largest SST 
anomalies resulted from sustained reversals of the wind 
direction. 

The results of a single synoptic calculation dur- 
ing the cooling season were presented in Section 6. This 
calculation showed that the inclusion of mean and anomalous 
advection resulted in an improvement of the accuracy of the 
model’s SST predictions. However, for the case studied, 
mean fluxes were used to drive the model and the anomalous 
currents v/ere computed using observed rather than predicted 
winds. In order to determine the model’s usefulness for 
predicting real time SST variations, predicted fluxes and 
winds must be used, and therefore the model should be coupled 
to a predictive atmospheric model for this purpose. 


The study has demonstrated the importance of 
horizontal advection in an ocean mixed-layer model. In 
coupling a mixed-layer ocean model to an atmospheric gen- 
eral circulation or numerical weather prediction model, it 
is clearly desirable that horizontal advection he included 
in the ocean as well as in the atmosphere, if the coupled 
model is to he useful for extended and long-range predic- 
tion, It is also apparent, however, that over a period of 
the order of one month or more, errors in the predicted sur- 
face wind fields may generate errors in the SST fields, 
which could result in a positive feedback of error into the 
atmospheric prediction. Planned experiments with a coupled 
ocean-atmosphere model should help to determine how serious 
this problem is likely to be with presently available models. 
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